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1 Introduction 

In recent years an abundance of knowledge on the correlation functions of integrable 
lattice systems has been obtained. Most of the development took place for the static 
correlations of systems with spin- 1/2 representations of the algebra su{2) including 
various anisotropic and inhomogeneous deformations. 

The developments started with the discovery of multiple integral representations for 
the density matrix of a finite chain segment [15, 19,20,27] and the observation [8] of the 
factorization of these integrals into sums over products of single integrals. With [4, 10] 
it became apparent that existence and factorization of multiple integral formulas are 
not tied to special properties of ground-states of infinite systems that often have higher 
symmetries. This holds even for finite temperature in the thermodynamic limit as well 
as for the ground states of finite chains motivating the research [3, 6, 7, 21]. Multiple 
integral formulations led to ab initio calculations of the asymptotics of ground state 
correlators of the XXZ chain [26]. Results for the closely related scalar Bose gas and 
for the Sine-Gordon model were obtained in [22,25,32,33]. 

Meanwhile the 'hidden Grassmann structure' was identified [6, 7] and made it 
possible to prove the complete factorization of general static correlation functions in 
terms of elementary nearest-neighbour functions under very general conditions [3,21]. 

In our view the most interesting question is if and how the investigations can 
be extended to models based on representations of higher rank algebras or to those 
based on higher level representations of su(2). We believe that the latter goal is 
in closer reach, but still very challenging. In particular, on the basis of [21] we 
expect algebraically transparent and physically relevant expressions for higher-spin 
generalizations constructed by means of the fusion procedure [30, 31]. 

Some older results of multiple integral type exist for the ground state and were 
derived in the ^-vertex operator approach in [9, 18] and in the combinatorial application 
of the algebraic Bethe ansatz in [24] for the isotropic case which was generalized 
to the zero field XXZ-case in [11]. The generalization of [24] to finite temperature 
and field was undertaken in [16]. Unfortunately, the resulting formulae are rather 
intricate for applications. In the present publication we undertake an effort to directly 
construct algebraic expressions for the finite temperature density matrices containing 
all information on correlation functions restricted to finite segments of infinitely long 
chains. Finite temperatures are treated by suitable 'lattice path integral' formulations, 
i.e. mapping the Id quantum chain to a suitable 2d classical system. Physical quantities 
are obtained by a combination of algebraic techniques (T- and ^-systems of transfer 
matrices) and analytical means (functional equations). 

A basic notion in the 'lattice path integral' approach to finite temperatures is the 
quantum transfer matrix (QTM) [35] acting in chain direction and being the 'quantum 
analogue' of quantum chains to the 'classical transfer matrix' of the classical Ising spin 
chain. The QTM satisfies the so-called ^-system from which by algebraic and analytic 
means the well-known thermodynamical Bethe ansatz (TBA) equations can be derived 
[23, 28]. Another way of analysis makes use of a finite number of auxiliary functions 
which satisfy a closed set of functional equations. This approach was developed for 
the spin- 1/2 case of XXX, XXZ and XYZ [29] and is probably the least canonical 



formulation, but highly efficient especially in applications. The generalization to the 
higher spin case was not clear before the work [34]. As we shall see below the auxiliary 
functions introduced in [34] are indeed most useful also for the explicit computation of 
correlation function of higher spin systems. 

In this work we derive algebraic expressions for the correlation functions of the 
integrable su{2) invariant spin-1 chain at arbitrary temperature, but zero field. We will 
be explicit for the 2-site case and give results for 3 sites. As the reader will see, the 
generalization to more sites and higher spin systems is viable. 

The paper is organized as follows. In section 2 we recall the construction of the 
Hamiltonian, the fusion procedure and fundamentals of thermodynamics. In section 3 
we introduce auxiliary density matrices and computational 'tricks', especially functional 
equations and the algebraic Bethe ansatz for a QTM with modifications yielding a 
'generating function' for the fundamental nearest-neighbour correlator in terms of which 
any other static correlation function factorizes. The fundamental nearest-neighbour 
correlator is characterized in section 4 in terms of useful integral equations. Explicit 
results in the low- and high-temperature limits are given. In section 5 we present 
concrete results for 2 sites. Section 6 is devoted to the complete computation of the 3 
site density matrix. A conclusion is given which is followed by appendices containing 
details of our construction. 

2 Fundamental thermodynamical and integrable structures 

2.1 Hamiltonian 

The Hamiltonian of the integrable isotropic spin- 1 chain on a lattice of L sites is 

/ ^ 

H=-Y,[Sn-i-Sn-{Sn~l-Snf], (1) 

^ n=l 

with for instance periodic boundary conditions Sl+i = Si. The spin components 5" act 
locally as standard spin-1 operators, and antiferromagnetic exchange, / > 0, is assumed 
throughout the paper. 

The Hamiltonian (1) was first obtained in a more general anisotropic form in [39]. 
Shortly later it was constructed by means of the fusion procedure [30, 31]. The ground 
state and the elementary excitations were studied in [38], and an algebraic Bethe ansatz 
and the thermodynamics within the TBA approach were obtained in [2]. 

2.2 Integrable structure 

The model can be constructed by means of the fusion procedure [30], starting from the 
fundamental spin-^ /?-matrix 



/?li''l(X) 



/I \ 

b{X) c{X) 

c{X) b{X) 



^(^)-I^' ^(^)-xf2i' ^^^ 



which we think of as an element of End(C^ (X" C^) (superscripts 1 indicate the spin- 1/2 
aka level 1 representations). It satisfies the Yang-Baxter equation 

As usual the /?L' in this equation act on the jth and kth factor of the triple tensor 

product C"^ (g)C'^ (g) C^ as R'^^'^' and on the remaining factor trivially. R^^'^' is normalized 
in such a way that 

/?[i.i](0)=pW, (4) 

where P^^' is the transposition of the two factors in C^ (S> C^. We say that /?['''! is regular. 
At the same time ^[^'^1 = p[i]/?[i'il satisfies the unitarity condition 

^[i,i](X)^[i-il(-X) = l4, (5) 

with 1„ denoting the nxn unit matrix. 

A further property of R^^-^', which is at the heart of the fusion procedure, is its 
degeneracy at two special points, 

where P^ are the orthogonal projectors onto the singlet and triplet subspaces V^"^ ,V^'' C 
C^ (g) C^. Due to (3) and (6) we have the important relation 

P23/?i'/] (X)/?S'/] {X + 2i)P+ = , (7) 

meaning that P^g^^ (?t)PiV' C^ + 2i) leaves C^ V (') invariant. 

Be 5 : C^ C^ — )■ C^ the projector onto the triplet space we write the fused R- 
matrices as 

/?[i'2l(X) = S23R^n\X)R^f{X + 2i)S'2, , (8a) 

/?P>i] (X) = SnRu^^ {X - 2i)43^' {X)S[2 , (8b) 

/?P.2](X) = SnS34R['4'\^-2i)R^i,'\X)R^^,'\X)R^^,'\X + 2i)SUS[2 (8c) 

acting on C^ (X" C^, C'^ (gi C^, or C^ (^ C^, respectively. Superscrips 1 and 2 refer to 
level 1 and level 2 (spin-1/2 and spin-1) representations. Combining the Yang-Baxter 
equation (3) and equations (7) it is easy to see that 

Rf^' •'■^■^1 {X - iu)Rf^' '^''^ (X)4^^^''''^l iiu) = Pg^''-'^1 {iu)Rf^' '"-^'l {X)Rf^' -^''^ (X-iu), (9) 

where i'y = ^ , 1 for j = 1 , 2, 3. 

In particular, /?P'^] is a solution of the Yang-Baxter equation. With P'-'^' denoting the 
transposition on C^ (S> C^ and pP'^1 = /'P]^[2,2] ^g have the standard initial condition 
and unitarity in the following form 

F'21(X)F'21(-X) = l9. (10b) 



It follows with (10a) that /?P'^1 generates the Hamiltonian (1), 

n=l 

2.3 Physical density matrix 

In [14] we have set up a formalism which enables us to calculate thermal correlation 
functions in integrable models with /?-matrices fulfilling (10a). It is based on the so- 
called quantum transfer matrix [35] and its associated monodromy matrix which are 
directly related to the statistical operator. 

The Hamiltonian (1) preserves the total spin 5" = Y}j=\S^. For instance, the 
magnetization in z-direction is a thermodynamic quantity, and the statistical operator 

Pi(r,/j) = e-P^-2«^ (12) 

describes the spin chain (1) in thermal equilibrium at temperature T{= 1/P) and 
magnetic field h. 

The free energy per lattice site /(r, h) determines the macroscopic thermodynamics 
of the model. For the general integrable spin-^ Heisenberg chain the free energy was 
calculated in [34] . 

The free energy is related to the partition function of the statistical operator Z := 
trp^. The normalized statistical operator 

Dl := ^PL, (13) 

is known as the density operator of the total system. 

It is convenient to define the reduced density matrix of a finite chain segment [l,m] 
by tracing out all other degrees of freedom 

The reduced density matrix is well-defined even in the thermodynamic limit L — ;• oo. 
With D[i^] we can calculate the expectation value of any local operator that acts trivially 
outside the finite segment [l,m]. 

For any integrable model, whose /^-matrix does not only satisfy the Yang-Baxter 
equation, but also the regularity and unitarity conditions (10), we can approximate the 
statistical operator pi{T,h) of the L-site Hamiltonian using the monodromy matrix of 
an appropriately defined vertex model with L vertical lines (1, . . . ,L) and A^ alternating 
horizontal lines {\,. . . ,N with N even). This fact was exploited many times in the 
calculation of the bulk thermodynamic properties of integrable quantum chains, in 
particular, in case of the higher-spin integrable Heisenberg chains [34]. In [14] it was 
noticed that the same formalism is also useful for the calculation of thermal correlation 
functions. Following the general prescription in [14] we define 



Tf^iX) = ^I^I^§(K + iu)I^iiu-X)... 



...l^^(X + iu)l^f''(iu-X), (15) 



where u := —J^/N and ti indicates transposition with respect to the first space in a 
tensor product. This monodromy matrix is constructed in such a way that (see [14]) 



2 ^ . . .. ,x /I 



^I...n{t!'^ (0) . . . ri'l (0) } = 1 - ^ I {Jhn-i,n - 2hSl) + 



n—i 



(16) 



Hence, the statistical operator pi can be approximated by pNi{T,h) := 
trj jv|ri'^'(0) . ..t[^\0)\. In the so-called Trotter hmit of A^ ^ oo we have Pn,l -^ Pl- 
The transfer matrix 

^Pl(X)=tr,-7}Pl(X) (17) 

is commonly called the quantum transfer matrix. We shall recall below how it can be 
diagonalized by means of the algebraic Bethe ansatz [34]. Quite generally it has the 
remarkable property that the eigenvalue A^l (0) of largest modulus of ?P1 (0) (we call it 
the dominant eigenvalue) is real and non-degenerate and is separated by the rest of the 
spectrum by a gap [35-37] . It can further be shown that 

f{T,h) = -T lim InAPl(O) . (18) 

Thus, the dominant eigenvalue alone determines the bulk thermodynamic properties of 
the spin chain. 

Owing to the fact that /?P'^1 satisfies the Yang-Baxter equation the transfer matrices 
f PI (X) form a commuting family, 

[?P](X),fPl(A/)]=0. (19) 

It follows that the eigenvectors of t'-^' (k) do not depend on X. Let |*P) denote a nor- 
malized eigenvector belonging to the dominant eigenvalue AP1(0). We shall call it 
the dominant eigenvector. It is unique up to normalization and is an eigenvector of 
?Pl(/'i) with eigenvalue AP1(X) = (*F|?Pl(/'i)|*F). In [14] it was pointed out that such an 
eigenvector determines all static correlation functions at temperature T and magnetic 
field h. In particular, it determines the density matrix (14) of any finite segment [l,m], 

n (Th^ i- (>y|rPi(o)^---^rP](o)|>y) 

DuM (T, h) = hrn^ (^[21(0))"' • ^20) 

For technical reasons it is better to consider a slightly more general expression 
than the one under the limit, by allowing for mutually distinct spectral parameters ^y, 
j = \,...,m, instead of zero. Setting ^ = (^i , . . . , ^m) we define 

the inhomogeneous density matrix at finite Trotter number. Then 

D[i,„](r,/i)= lim lim dPI(^). (22) 



3 Modified lattice model and auxiliary density matrices 

The physically interesting correlation functions are spin-1 correlators, however, the 
most fundamental objects are spin- 1/2 correlators resp. spin- 1/2 density matrices 

Al'l(i;,)...Al'l(i;„) 

Both types of density matrices are properly normalized, i.e. the total traces evaluate to 
1, because of the normalization terms in the denominators. Taking traces oi D^J^ turns 
the monodromy matrices in the numerators T^-''{^) into transfer matrices r-''{^) which 
evaluate to the eigenvalue AI^T (^) when acting on the eigenstate |*P). 

The density matrix D^l (^i , . . . , ^m) is obtained from dI^I (^i , . . . , ^„) by fusion in the 
case n = 2m and symmetrizations over neighbour pairs of quantum spaces carrying 
spectral parameters ^2/t-i =^k — h ^2k = ^yt + i- This procedure yields the unnormalized 
DPl(^i,...,^m)- The normalization factor is given by the ratio of denominators occuring 
in (23), (21) 

- AWfo-i)AWfo + i) 
l\ AP1(^,) • ^"^^ 

In this paper we follow the strategy of (discrete) functional equations for finite tempera- 
ture [1] to calculate the spin-1 correlators. This approach is actually formulated and 
viable for finite Trotter number A'^. It yields A'^ many equations for the density matrix 
dP1(^i, ...,^m) as function of one of its spectral parameter arguments, let us say the 
last one ^m- Together with the asymptotic condition lim^^^^D'-^'{^i, ...,^m-b^m) = 
DPl(^i,...,^m-i)'Xil (for zero field) it fixes the dependence of dP1(^i, ...,^„,) on ^^ 
if the data dP1(^i, ...,^m-i) are given and suitable analyticity conditions hold. The 
analyticity conditions in the spin-1/2 case [1] are simple: The entries of the density 
matrix for a system with Trotter number A^ are multivariate polynomials of degree A^ 
divided by the known multivariate polynomial A(^i)...A(^,„). Hence, the A'^ + 1 many 
equations (linear independence was proven in [1]) fix the A'^+ 1 many coefficients. 

For our case at hand two different approaches are in principle conceivable: (i) the 
spin- 1 object D^l (^i , . . . , ^m) is tackled directly, or (ii) the spin- 1/2 object d'^I (^i , . . . , ^„) 
is treated first and then fusion is applied. Unfortunately, both approaches present 
problems: (i) the object d[^](^i,...,^,„) consists of multivariate polynomials of degree 
2A^ (!); (ii) for D[^1(^i,...,^„) which consists of multivariate polynomials of degree A^ 
the functional equations can not be derived. In both cases though, the problems can be 
overcome. 

3. 1 Functional properties of D I^l ( ^ i , . . . , ^^ ) 

In this paragraph we review the basic constructions for setting up (discrete) functional 
equations for reduced density matrices of integrable vertex models on semi-infinite 
(oo X A^) square lattices. It is possible, and for later purposes convenient, to start 
with a system with A'^ many different, but fixed spectral parameters uj on the A'^ many 
horizontal lines. The (unnormalized) density matrix on a finite sequence of sites has 



matrix elements that are partition functions of 'sliced lattices'. For notational simplicity 
we write D(^i,...,^m) for DPl(^i,...,^m) wherever possible. The local /^-matrix enjoys 
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Figure 1: Graphical illustration of the matrix element D^'j'''.a" (^i)---)^m) for m = 
4. The lattice is semi-infinite with N many infinitely long horizontal lines carrying 
independent spectral parameters u\ , ...,un. 



standard initial condition and crossing symmetry which allows to perform 'lattice 
surgery' as described in [1]. If one of the spectral parameters ^1, ...,^m> let us say the 
last one ^m, on the vertical lines with cut happens to be identical to any of the spectral 
parameters on the horizontal lines, the matrix D(^i , . . . , ^m) rnaps to D(^i , . . . , ^^ — 2i) 
upon a linear map A acting on the (2^'" -dimensional) space of reduced density matrices. 
This yields the (discrete) qKZ-type functional equation 
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a 6 6 ^4 - 2i 6 6 6 ^4 

Figure 2: Graphical illustration of the functional equation (25). 

D(^i,...,^^-2i)=A(^i,...,^„,)(^b-,^m). 



(25) 



The map A depends on the spectral parameters ^1, ...,^m and is given by a product of 
2m — 2 many /^-matrices. 

In addition we have the asymptotic condition 



lim D(^i,...,^^_i,^m) =D(^i,...,^„_i)®l. 



(26) 



The under-determinacy of the functional equations for the density matrix as function 
of a single variable is simply resolved by considering the full dependence on all 
arguments and invoking the intertwining (symmetry) relations 

Rk,k+^i^k,^k+l)D{^l,...,^k,^k+l,---,km)Rkk+l(^k,^k+l) =D{t,i,...,^k+h^k,---,km)- 

(27) 
Let us be explicit for the 2-site case which is non-trivial and will be analyzed in 
detail in the next sections. 

^1^2(^1, ^2)D[2l(^i,^2)^-j(^i, ^2) =dP1(^2,^i). (28) 

Here we explicitly remind of the spin-1 (level 2) case by the superscript 2. Due to su{2) 
invariance some simplification occurs which is not essential for the following reasoning, 
but leads to some notational simplification. The 2-site density matrix as well as the 
intertwiner ^12 are superpositions of projectors onto singlet, triplet and quintuplet 
spaces. All of these operators commute. Hence, ^12 and ^j^2 o" the left hand side 
of (28) drop out and we obtain D^l (^1 , ^2) = D^^'^ (^2,^1)- Each of the entries of this 
density matrix is a symmetric polynomial of two arguments of degree 2N (divided by 
the product of eigenvalues A^l (^i)API (^2))- Each of these polynomials has {2N + 1)^ 
coeffcients. Due to symmetry, only {2N+ 1)(A^+ 1) coefficients are independent. There 
are equations for D^l (^1 , ^2) for arbitrary first argument ^1 and for A^ + 1 many discrete 
values of the second argument ^2 = m 1 ,-••,"«, °° which amount to {2N + 1 ) (A'^ + 1 ) many 
equations for the coefficients. (We have checked linear independence for characteristic 
cases.) 

The above reasoning shows that the solution to the functional equation (25), asymp- 
totics (26) and symmetry (27) is unique for the analytic condition that each matrix 
element is of the type 

PK'jl) •••iSmJ .»q. 

AP](^0-AP1(^^) 
where pis a multivariate polynomial of degree 2N. Finding the solution is a different 
story and is based on a suitable ansatz. The result will be given in the next sections. 

The above reasoning was done for pairwise different spectral parameters on the 
horizontal lines of the lattice. This allowed for a simple counting of independent 
equations. In the case of degenerate values, the equations can be formulated for the 
functions and for certain derivatives yielding the same counting. In the following we 
are dealing with systems with N/2 many spectral parameters iu and N/2 many spectral 
parameters 2i — iu. Clearly, the choice of these parameters can be relaxed to have 
pairwise different numbers that still approximate the statistical operator. 

3.2 Functional properties of D^^^ (^1 , ..., ^«) 

From a fundamental point of view it would be highly desirable to tackle the density 
matrix D['](^i,...,^„) (vertical spin- 1/2 lines embedded in spin-1 lines). There are no 
clear functional properties like in the D^l case as here the mixed /^-matrices do not enjoy 
the standard initial condition. However, D[^1(^i, ...,^„) may be obtained in a certain 
limit of a system with (clearly) n many vertical spin- 1/2 lines, 2N many horizontal 
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^1 ^2 6 ^4 "0" 

Figure 3: Graphical illustration of the auxiliary lattice. If £ = is chosen and sym- 
metrizers at neighbouring pairs of horizontal lines with spectral parameters Uj it i are 
attached, the lattice will turn into a realization of D^^l (^ 1 , . . . , ^„ ) . In the main text, an 
argument is given to as why the same result appears in a simple £ — )■ limit without 
symmetrizers. This argument is based on the infinitely many column-to-column transfer 
matrices with spin-1 auxiliary space, cf. the boxed object in the figure. 



spin- 1/2 lines and infinitely many vertical spin-1 lines, see Fig. 3 for the case n = 4. The 
naive idea is to use a direct fusion construction based on a lattice like the one shown in 
Fig. 3 with £ = and symmetrizers applied to the horizontal spaces. Next one would try 
to write down 2N many functional equations of the type (25) as now the /?-matrices of 
the vertical spin-1/2 lines and the 2N many horizontal spin-1/2 lines enjoy the standard 
initial condition and crossing symmetry. Unfortunately, the spectral parameters on 
the horizontal lines are not independent. Only A^ many meaningful equations can be 
derived, the other N many attempts result in uncontrolled 0/0 expressions. 

The system shown in Fig. 3 is well-defined and has a regular £ — )• limit: the 
column-to-column transfer matrices with auxiliary spin-1 space take their leading state 
close to the space of states that are symmetric in neighbouring pairs. In the limit £ — ;■ 
the leading state converges to a state in the pure spin-1 subspace (C^) of the total 

space (C^) . This is clear from the root pattern of the Bethe ansatz solution to the 
system Fig. 3. 

The benefit of dealing with the '£-regularized' system is the possibility to now 
formulate 2N many functional equations. These equations are valid for any small, but 
finite £ for the semi-infinite system. This is the program of the pure spin-1/2 case 
that leads to algebraic expressions of the density matrix in terms of nearest-neighbour 
correlators [1,21]. The algebraic structure is completely independent of £! Just the 
nearest-neighbour correlators depend on £, but have a regular limit. We therefore see 
that the factorized expressions also hold for the object d'^I (^1 , . . . , ^„) which is of interest 
to us. Now the computational strategy is to write down the entries of D['](^i,...,^„) 
for n = 2m in terms of the fundamental nearest neighbour function ^(^1,^2) that is 
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calculated in the next section. Then the symmetrization in the quantum spaces (vertical 
lines) is to be performed for ^2yt-i =^k — i,^2k = ^k + i{k=i,---,m). 

3.3 Bethe Ansatz solution 

For the calculation of the free energy (18) and the inhomogeneous density matrix (21) 
we need to know in first place the dominant eigenvector |*F) and the corresponding 
transfer matrix eigenvalue aP1(X). They can be obtained by means of the standard 
algebraic Bethe ansatz for the spin-^ generalized model (see e.g. chapter 12.1.6 of [12]), 
since, by the general reasoning of the fusion procedure [31], the quantum transfer 
matrix r"^' (k) can be expressed in terms of a transfer matrix with spin-| auxiliary space. 
For our purposes we consider a staggered monodromy matrix with spin- ^ auxiliary 
space and N many spin-1 quantum spaces and two spin- 1/2 spaces by 



ri^l (X + i) = eP'-^/?i;;;^2 (^ - A/ + m^N+t (a- + 5 - X - i) • 



Then, interpreting this monodromy matrix as a 2 x 2 matrix in the auxiliary space a, we 
define the transfer matrix and the quantum determinant (spin-0 fusion) as 

?W(X) = trrW(X), detqrW(X) = ^(rW(X-i)®rW(X+i))^', (3i) 

where U is the projector onto the antisymmetric state in the tensor product of the two 
auxihary spaces. It follows from (8) that 

rP](X) = 5(rW(?i-i)^rW(X+i))5', (32) 

where S is the projector onto symmetric states in the tensor product of the two auxiliary 
spaces. Taking the trace and using (31) we conclude that 

?[2l(X)=?[il(X-i)?W(X + i)-detqrW(X), (33) 

Hence, since detqr['l(X) commutes with T^^\X) [31], every eigenstate oi t^^\X) is an 
eigenstate of r^l (k) as well. 

The algebraic Bethe ansatz is based on the Yang-Baxter algebra relations 

^[^•ii(X-A/)(rW(X)0rW(A/)) = (rW(A/)^rW(X))^[''ii(X-A') (34) 

which follow from (9) and (30). Representing jl^l (k) by the 2x2 matrix 

^ ' \C{X) D{X)j ' 
the quantum determinant is 

detq rl'l {X) =D{l-i)A{l + i)-B{'k- i)C(X + i) . (36) 
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Defining the pseudo vacuum 

|0)=[(o)^(o)]''^^(?)®a) (37) 

we deduce from (30) that 

C(X)|0)=0, A{1)\0) = a{l)\0) , D{X)\0) = d{X)\0) , (38) 

where the pseudo vacuum eigenvalues a(X) and d(k) are explicit complex valued 
functions. Using the notation 

<\>±iX) = iX±iuf/\ « = -^. (39) 

which proved to be useful in [34], and 

(p+(X) :=X-^, (p_(X) :=X-^-5, (40) 

which arise from the additional horizontal spin- 1/2 lines carrying the spectral parameters 
/7 — i and /u — i + b (at the vertex with transposition), we can express the vacuum 
eigenvalues as 

^^^ ^_(X-3i)(p_(X-2i)' ^^^' ^+(X + 3i)(p+(:)i + 2i)- ^ ' 

Given the Yang-Baxter algebra (34) and the pseudo vacuum eigenvalues (41) the 
eigenvectors and eigenvalues of r^'Ck) can be obtained from general considerations 
(see e.g. chapter 12.1.6 of [12]). The dominant eigenstate \'^) of ?Pl(?i), in particular, 
can be represented as 

|»F)=B(Xi)...B(X^+i)|0), (42) 

where the set of so-called Bethe roots {'kj}^J^^ is a specific solution of the Bethe ansatz 
equations 

For the given set of Bethe roots {?ty}^^/ we define the Q-function 

q{X)=Y{(k-Xj). (44) 

7=1 

Then the eigenvalue of ^^1 {X) corresponding to |*P) is 

As for the eigenvalue of f PI {X) we conclude with (33) and equation (36) above that 

APl(X)=AW(X-i)AW(X + i)-a(X + i)J(X-i) (46) 

This eigenvalue and the Bethe ansatz equations (43) are the main input for the calculation 
of the thermodynamics of the spin-1 chain. In order to perform the Trotter limit the 
eigenvalue must be represented by means of auxiliary functions satisfying a finite set of 
nonlinear integral equations. This was achieved in [34]. 
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4 Integral equations and basic functions 

In this section we consider the evaluation of the largest eigenvalue of the generalized 
quantum transfer matrix by means of nonlinear integral equations (NLIE). This gives 
us the opportunity to introduce certain auxiliary functions that are also relevant for the 
factorized algebraic expressions of the density matrix elements in the next section. Our 
starting point is the expression for the dominant eigenvalue A^l (X) of the quantum 
transfer matrix together with the Bethe ansatz solution (43)-(46). In [34] the problem 
was solved within the more general context of the fusion hierarchy, and NLIE for the 
integrable isotropic spin chains of arbitrary spin were obtained. 

For the calculation of the free energy for spin 1 we will be dealing with three 
coupled NLIE for three functions b, b and y. We show the equations below in (48) and 
present the derivation in appendix B. (Most of the relations of this section are valid for 
arbitrary magnetic field h, however the main applications later on concern the zero field 
case.) For finite Trotter number A^ the functions b, b and y can be expressed in terms 
of the 2-function (44) and the functions <^± introduced in (39) (see appendix A). This 
defines them as meromorphic functions on the entire complex plane, but is inappropriate 
for performing the Trotter limit. In the NLIE, on the other hand, the Trotter number 
appears only in the driving term and the Trotter limit is easily obtained. For a discussion 
of some of the subtleties related to the Trotter limit and the definition of useful auxiliary 
functions see [17]. 

The functions b(X),b(X) and j(?i) are defined in (A.5), (A.7) and will be of particu- 
lar use in the neighbourhood of the real axis. For convenience we introduce the shifted 
functions 

be(X) = b(X-i£), be(X) = b(?i + i£), (47) 

and similar capital functions, !Be(X) = 1 + be(^) etc. Then the desired NLIE read 

'logy{X)\ /a,(X)\ (logY{X)\ 

logbe(X) = Ah{X) \+%*\ log^e(^) , (48) 

^logbe(X)/ \A^{X)) \log^e{l)J 

where (JC*g)i denotes the matrix convolution ^y/!°oodX'3C,j(X — X')gy(X'), and 

tanh(f(X-A' + i)) 



Ay{X) =l0g 



tanh(^(X-A(-5 + i)) 



cosh I (A — /^j 



AhiX) = -^+diu,X-ie), A^iX) = ^+diu,X + ie), (49b) 

d{uM = ^j Jk. '^^^^^ -. -^^. (49c) 

The integration constants {±h/T) are fixed by comparing the asymptotic values of both 
sides of (48) for |X| — )■ oo. The kernel matrix is given by 

Xil + ie) Xil-iE) 

%{'k)= {XCk-ie) J(X) -J(X + 2i(l-£)) I , (50) 

.X{X + ie) -J(X-2i(l-£)) J(X) 
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where 



%{X) 



, n^) = 



dk e 



-\k\-ikX 



'-co 271 2 cosh ^ 



(51) 



4cosh7lX/2 

and 3" can be expressed in terms of the digamma function *?, see (56). 

The derivation of integral expressions for the eigenvalues is involved, we defer the 
details to appendix B. The results are 



lOgAW(X): 



r A^e-31^1 cosh ((m + \)k) ^_ 

(XK ■ '^ 



2\k\co&\\{k) 



■e-'*^^-log[^+(X + 3i)^_(X-3i)]+cst. 



+ 



-/ JV3^(X-V)|logF(V) + log 
respectively 

logAP](X)=log 



(p+(X-i)(p_(X + i) 
(p_(X-i)cp+(X + i) 



(52a) 



(^+ (X-2i)(^_ (X + 2i) (p+ (X-i)cp_ (X + i) 
(^_(X-2i)(|)+(X + 2i)(p_(X-i)cp+(X + i) 



+ log3;(X) (52b) 



Next we are interested in ^logA^^l for which we derive an integral expression in 
terms of functions satisfying a set of linear integral equations. The basic functions 
are the derivatives of the functions logy,log*Be,log!BE with respect to 8. We have the 
identities 



Gy{X,^i):=^\ogY{X) 



35 



8=0 



35 log)' W 



8=0 



[l+y-\X)]\^=oGy{K^i), 



(53a) 



Gi,(X,A/):= ^logSe(^) 



%(?i,A/):=^log«e(X) 



8=0 



8=0 



36 

3 
35 



log be (k) 
log be (X) 



= [l + b-HX)]|8=oGb(^,A'), 
8=0 

(53b) 

[l + b;'(X)]|8=oG|^(X,A(). 

(53c) 



8=0 



These satisfy the set of linear integral equations 



[\+y-\X)]Gy{K^i)\ f-i 
[l + K\X)]G^iX,^)) \ 



^coshf(X-;u) 






^Gy{X,fi)\ 

+ X*\ Gb(X,ij) 
\Gk{X,^)J 



(54) 



Next we obtain an explicit expression of the eigenvalue's derivative in terms of the G 
functions 



35 



logA[il(X) 



= f dX' 7—- — 7r-:vlGy(X',^) + — ^.5 I (55a) 

^0 i-cc 4cosh(|(?i-V)) I ^ (X'-^y + lj 

r°° 1 
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for I Im(?i)| < 1. Next we note a useful identity for 3^ defined in (51) 



n^) 



1 



"¥ 



:l 



+ ^{-X 



871 1 V 4 
where *P is the standard digamma function 



"¥ 




'i' x) 



dx 



iog(r(x)) , r(x) 



f-^e-'dt 



We remind of the functional equations 



and 



»F(l-x)=»F(x)+7Ictg(7lx) , x/ 



»P(x+l)=»F(x) + 



1 

2' 



The eigenvalue's derivative satisfies the following 2-point equation 



36 



lnAW(X + i) 



+ ^lnAW(X-i] 
5=0 ^^ 



Gy{X,^) + 



1 



1 



8=0 



X — fd — i X — /u + i 



(56) 
(57) 

(58) 
(59) 

(60) 



The right hand side does not simplify in any obvious way to elementary expressions due 
to the occurrence of Gy{X,^). In the next subsection we will derive a simple 3-point 
equation involving just derivatives of A^^]. 

Due to su{2) invariance the density operator D[^1(X,/^) for two neighbouring spin- 
1/2 spaces with spectral parameters X,^ has the following representation 



dW(X,a<) 



^-^(0(X,^)]1 + ^(0(X,A<)PW 



(61) 



where P^^^ is the permutation operator of neighbouring spin- 1/2 objects and a)(X,;u) is 
some (symmetric) function. In Appendix C we show that 

For later purposes, it is convenient to use another closely related function 

2i 



(62) 



a{X,n) :-- 



2i 



(X-A^)2+4) 



tr{DW(X,/i)pW} 



38 



InM 



Hk^)} 



+ 



2i 



8=0 (X-/j)2+4 



(63) 



(X-/i)2+4 

An interesting and useful observation concerns the trace of P'^^' over the symmetric 
subspace for arguments (X,/^) chosen as (X — i,X + i) which according to the fusion 
principles yields 

A2(X) 1 



3 1 

- + -oi{X-i,X + i) 



(64) 



Ai(X + i)Ai(X-i) l+y{l)-^' 

We have the explicit expression (62) which seems to yield j forco(X — i,X + i). However, 
^ InAl'l (X;/7) develops poles for X — fj = ±21 which can be evaluated on the basis of 
(55a) and (54) in full agreement with (64). 
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4.1 Functional equations for the basic functions 

From the definition (A. la), (A. lb) of the functions Ai and A2 we find the relations 

Ai {I - i) Ai (X + i) = ^(X - 3i)^(X + 3i)cp(X - 2i)cp(X + 2i) + A2 (l) 

= (^{l - 3i)^(X + 3i)(p(X - 2i)(p(X + 2i) • ¥{1) (65) 

A2(X-i)A2(X + i) = ^(X-2i)^(X-4i)^(X + 4i)^(X + 2i)x 

X cp(X - i)(p(X - 3i)(p(X + 3i)cp(X + i) +^{X)A(X) . (66) 

We divide the first equation by (p(X — 2i)(p(X + 2i) and the second equation by a product 
of that function with X replaced by X it i 



cp(X-2i)cp(X + 2i) ^^ '^^ ' (p(X-2i)(p(X + 2i) 

(67) 

A2(X-i) A2(X + i) 



^(X-2i)^(X-4i)^(X + 4i)^(X + 2i)+ 



(p(?i-3i)cp(X + i) cp(X-i)cp(X + 3i) 

, m Am 

cp(X-i)(p(X-3i)cp(X + 3i)(p(X + i) ^'' 

(68) 

From the last equation we see that 



as (p(X-3i)(p(X + i) 35 (p(X-i)(p(X + 3i) ^^ ^ 

since for V% with (|)(X) = the right hand side of (68) is absolutely independent of 5. 

a A2(?i) _(p(X-2i)cp(X + 2i) a A2(X) 



35 (p(X-2i)(p(X + 2i) A2(X) a5cp(X-2i)(p(X + 2i) 

_ cp(X-2i)cp(X + 2i) a Ai(X-i)Ai(X + i) 
~ A2(X) a5(p(X-2i)cp(X + 2i) 

^ Ai(X-i)Ai(X + i) 3 Ai(X-i)Ai(X + i) 

A2(X) 35 ''(p(X-2i)(p(X + 2i)' ^ ^ 

where for the second identity we have used (67) where the first summand on the right 
hand side is independent of 5. Next, inside the logarithm, we replace Ai by At^l and the 
logarithmic derivative by the function Q. (63). We then find for the right hand side of 
(70) 

Ai(X-i)Ai(>i + i) /^ 3i i on ; ^ on 1 ; ^ 



A2(X) V(^-A')' + 9 (X-A')' + l 

(71) 

The ratio of A functions can be simplified by use of (64). With 
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we obtain 

Finally, the functional equation (69) takes the form 

Q.(X-2i,^) +Q{X,fi) - o(X- ^-i) QiX,^) +D.{X + 2i,fi) - o{X- ^ + i) _ 

N{X-i) ^ N{X + i) ~ ' 

(74) 
for all X being zeros of (|)(X). 

4.2 Low and high temperature Umits 

Low temperature limit, zero field 

From (48) we read off for zero field h = the zero-temperature limit T — ;■ of the 
auxiliary functions y, be and bg 

lim v|s-n = 1 ) lim bgU-n = = lim bgU „ • (75) 

And from (54) we obtain in the same limit 

Hence the zero-temperature limit of the functions Q. and co (63) is calculated directly 
from (55b) yielding 

lim a)(V) = \- \{iX-^f + A) \U> (-^(X-A')) +'i' (^(^-A'; 

///g/i temperature expansion, zero field 

We expand the auxiliary functions in powers of \/T 



3'Wl8=o=3'o+£(^) yk{X), (78a) 

be(^)l8=o = t'o + f:(^) bkiX), (78b) 

^eWl8=o = ^0 + £ (^) b,(X) , (78c) 

and solve the NLIE iteratively. The O"^ order is solved by the constants jo = 3, bo = 
b() = 2. The NLIE for the auxiliary functions up to P' order linearize in y^ {X), bi (X), 
bi (X) and can be solved by Fourier transformation. This can be repeated over and over 
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again: the NLIE linearize in the next (unknown) order where the (known) lower order 
terms appear in non-linear combinations as driving terms. We restrict ourselves to the 
first order calculation 

n^\ - 128 7 



(4+X2)(i6+X2)r 



+ o(i/n, 



^e(^)l5=o 



32i(4 + i(X-i£) + (X-i£) 



/ 



b^{X) 



8=0 



:2 + 



((X-i£)2 + l)((X-i£)2+9)(X-i£ + 5i)r 

32i(4-i(X + i£) + (X + i£)2) / 

((?l + i£)2 + l)((X + i£)2+9)(X + i£-5i)r 



+ 0{\/T^ 



+ 0{\/T^ 



(79a) 
(79b) 
(79c) 



Using this strategy also for the auxiliary functions Gy, G^, G^, the high temperature 
expansion for the logarithmic derivative of the eigenvalue with respect to 8 leads to 



log All] (X) 



38 
Thus we have 



s=o 



(0(X,A/) 




+ o(i/r2 



•+X2)(9+A(2)r 



(80) 



(81) 



5 Computation of the 2-site density matrix 

Due to rotational invariance the density operator of two neighbouring spin-1 sites has 
to take the form 

dP1(^i,^2) = Pi(^i,^2)1+P2(^i,^2)/'P1 + P3(^i,^2)/'("\ (82) 

where pP] is the permutation operator and p(") := (3 x) the projector onto the singlet. 
Three coefficients pk,k= 1,2,3, have to be determined. 

The general functional equation of qKZ type specialized to the 2-site case and 
written in terms of the above introduced coefficients reads 




^M^l-^2) 



Pl(^l,^2) 

P2(^l,^2) 

,P3(^1,^2) 



(83) 



Using the shorthand notation ^ := ^i 

' (5-2i)( 



(5-2i)(5+2i)(5-4i)(^+4i) 



^2 the matrix L(^) is given by 

4i5(4-4i5H^) 4i^ \ 

Mi) \ 



(^-2i)(5+2i)(5-4i)(^+4i) 

-4i^(8-6i^H^) 

(^-2i) (^+2i) (^-4i) (5+4i) (^-2i) (5+2i) (5-4i) (^+4i) 
I 16(12+^^) 64-16i^-4^^-8i^^+^^ 

\ (5-2i)(^+2i)(^-4i)(^+4i) (^-2i)(5+2i)(5-4i)(^+4i) 



(5+2i)(5^ 
^(^-2i) 
(5+2i)(^+4i 
-4i(^-2i) 
(5+2i)(^+4i)/ 



(84) 



The functional equations for the three coefficients may be disentangled to some degree 
by the transformation 



(^1,^2 



16 



^2+4 



i^'- 



w- 




(85) 
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In fact, F is fixed by tlie trace condition of the density operator 

F(^i,^2)=9pi(^i,^2) + 3p2(^i,^2)+3p3(^i,^2)=trDPl(^i,^2)^l. (86) 
Hence we find tlie simpler looking functional equations 

l\ / 1 \ (l 



shift 



\- 



1 



256;(^-i) 



3(^+2i)(^-2i)2(^+4i) 





(^-2i)(^+4i) 
^(^-6i)(^^-2i^-4) ^^(^-6i)(^-4i) 

■ (^-2i)2(^+2i)(^+4i) (^-2i)2(5+2i)(5+4i) / 







(87) 

where the functions G,H on the left hand side are taken at shifted arguments (^i — 2i, ^2). 
and those on the right hand side at (^1 , ^2)- The inverse transformation is 



Pi(^i,^2r 

P2(^l,^2) 
,P3(^1,^2) 



/ 5^^+36 
' 45(^2+4) 
64 



V 



' 30(^2+4) 

3^2-20 

"45(^2+4) 60(^2+4) 

16 3^2^20 

45(^2+4) 60(^2+4) 



^^+6 \ 
15(^2+4) \ 

3^2+28 
30(^2+4) 

3^^+8 ■ 
30(^2+4) / 



G(^l,^2) 
,^(^1,^2) 



(88) 



Below, we show that the explicit and unique solution for the functions py(^i,^2), 
j= 1,2,3, is 



Pi (^1,^2 



1 



1 



iV(^i)A^fe) 180^2(4 + ^2 
3 



■4(5^^(4 + ^^)) 



-(5e(4 + ^2)-96)((0(^i-i,^i+i)(0fe-i,^2+i)) 

-3(4 + ^2)((^ + 2i)(^-6i)(0(^i-i,^2 + i) + (^-2i)(^ + 6i)(0(^i+i,^2-i)) 

+ 3^2 (16 + ^2) ((0(^^_i^^^_i) + (0(^j+i^^2+i)) 

+ (6 + e) (^2 (16 + ^2) (0(^1 -i,^2-i)(0(^l+i,^2 + i) 

(89a) 



(4 + e)'(0(^l-i,^2+i)(0(^l+i,^2-i))}, 



P2(^l,^2) 



1 



1 



iV(^i)A^(^2) 180^2(4 + ^2) 

12(5^2_28)((o(^j_i,^^+i) + (o(^2-i,^2 + i)) 

+ S{5^^-2S)i(0i^,-i,^,+iM^2-h^2 + i)) 

-3(4 + ^2)((^ + 2i)(^+14i)(0(^i-i,^2 + i) + (^-2i)(^-14i)(0(^i+i,^2-i)) 

-12^2 (16 + ^2) ((0(^1 -i,^2-i) + (0(^i+i,^2 + i)) 
1 



28 + 3^2 ^2 16 + ^2 «(^i_i, ^2- i)(o(^i+i, ^2 + i) 



P3(^l,^2) 



(4 + ^2)^^(^^_i^^^+i)^(^^+i^^^_i^ 
1 1 



(89b) 



A^(^i)A^(^2) 180^2(4 + ^2 
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-12 (8 + 5^2) ((0(^1 -i,^i+i) + (ofe-i,^2 + i)) 

-8 (8 + 5^2^ ((0(^1 -i,^i+i)(ofe-i,^2+i)) 

+ 12(4 + ^2) ((^-i)(^ + 2i)(0(^i-i,^2 + i) + (^-2i)(^ + i)(0(^i+i,^2-i)) 

1 



(89c) 



2(8 + 3^2)(^2(j6^^2)^(^^_i^^^_i)^(^^^i^^^^i) 

Note that the function A'^(^) appearing in the denominator is defined in (72). 
Proof: The symmetry of the function co 

(0(^1,^2) = (ofe,^i), (90) 

directly leads to the symmetry of the coefficients Pit(^i,^2)- The asymptotics of the 
function co 

lim (0(^1,^2) =0 for any 7 = 1,2, (91) 

as obtained from (26), for instance, entails the asymptotics of the functions Pit(^i , ^2) 

limp,(^i,^2) = (^ ^°'^ = ^' (92) 

5;^~ I else. 

The functions p/:(^i,^2) as solutions of the functional equations were constructed in 
the following way: (74) is independent of the second spectral parameter, hence also the 
function 



Q(^i+i,^2+i)+n(^i-i,^2 + i) + n(^i+i,^2-i) + ^(^i-i,^2-i) 

-o(^i-^2-i)-o(^i-^2 + i)} (93) 

satisfies (74), where C is some constant. Then it is easy to show that the function 

K^i,^2):=^^^^^^^^|^(^i+i,^2+i)^(^i-i,^2-i)-n(^i-i,^2+i)^(^i+i,^2-i) 



+a(^i-i,^2+i)o(^i-^2+i)+^(^i+i,^2-i)o(^i-^2-i)-o(^i-^2+i)o(^i-^2-i) 

(94) 

satisfies 

/j(^l,^2)-K^l-2i,^2)=o(^l-^2-ik(^l,^2) (95) 
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Define 

G(^i,^2) := ^(^^ + 16)(e + 4)g(^i,^2) , (96a) 

//(^i,^2):=^(e + 16)(^2^4)2M^i,^2) + ^(|) , (96b) 

then these functions satisfy (87). Applying the inverse transformation (88) the constant 
in (93) and (94) is fixed by the analyticity requirement for the functions pk{^i,i,2), 
leading to the cancellation of spurious poles at ^i = ^2 ■ Using the identity (64) or (72) 
respectively leads to the results (89). 

Finally, we consider the low and high temperature limits of Pfc(^i,^2) in the homo- 
geneous case ^1 =^2 = 0. With the results we easily calculate any nearest neighbour 
correlators like for instance the internal energy e = —T^^ if /T) by 



A=o = -^[2pi(0,0) + 3p3(0,0)]. (97) 



_3/ 

~y 

Low temperature limit, zero field 
By use of (77) we obtain 

12 19 14 14 

P,(0.0) = ---.^ p.(0,0) = -- + -.^ P,(0.0) = -- + -.^ (98) 

Seemingly no zeta function values appear in these expressions. At first sight, this looks 
qualitatively different from the case of the spin- 1/2 Heisenberg chain. Here, for the spin- 
1 chain, we find a 71^ expression instead of zeta functions at odd integers. However, this 

2 

can be rewritten in terms of Riemann's zeta function at the special argument ^(2) = ^. 
We will see below that this pattern continues and for the 3-site density matrix we obtain 
rational degree 3 polynomials in 71^ or equivalently linear expressions in ^(2), ^(4) and 
^(6) with rational coefficients. 
For the energy in the zero-temperature limit we simply get 

-7, (99) 

without any 71^ contributions. This result is in agreement with the literature, cf. [38]. 

High temperature asymptotics, zero field 
By use of (81) we obtain 

Pi(0,0):^^-^^ + 0(r-2), (100) 

P2(o,o)c^-^^ + o(r-2), (101) 

p3(0,0)::^^^ + O(r-2), (102) 



and for the energy 



U = -3-9^ + 0(r-^)- (103) 
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6 Computation of the 3-site density matrix 

Similar to the 2-site case in the preceding section, the density operator of a 3-site 
segment is presented in the form, 



15 

I 

a=l 



DPl(^l,^2,^3)=IPa(^l,^2,^3)/'a. 



(104) 



The projectors Pa (1 < oc < 15) are graphically represented in Fig. 4. The lines connect- 

12 3 4 5 

• • • 

• • • 





6 7 



9 
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10 
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12 



13 



14 



15 



Figure 4: 

ing upper and lower rows indicate Kronecker deltas. The lines connecting nodes a and 
b {\ <a <b <3)inthe same row represent matrix elements C{ia,ib)^ia+k,0 such that 



C{ia,ih) 



1 {ia,ih) = (0,0) 



(105) 



[-1 (/„/,) = (1,-1), (-1,1). 

This choice of projectors ensures the SU{2) symmetry of the density operator, 

GdP1(^i,^2,^3)G-i=dP1(^i,^2,^3), GgSU(2). (106) 

For later convenience, we also introduce renormalized coefficients 



P.(^l, ^2,^3) = n^(yP.(^l, ^2,^3) 



(107) 



where N{i,j) is defined in (72). 



23 



The result of our calculations takes a compact form in the homogeneous limit 

^1 = ^2 = ^3 = at r = (where p« = 8pJ: 



1879 3497 . 53 . 11296 . 

Pi = 71 H 71 71 

^ 432 1350 135 637875 

^ ^ 953 2 37 4 1043 1552 . 

po = p, = 71^ ^ 71^ -I 71° 

^^ ^^ 2700 675 2160 637875 
^ ^ 983 2 251 1592 ^ 1 . 

P4 = P^; = 71^ j 7t" 71^ 

^^ ^^ 900 144 212625 6 
197 2 374 4 13021 3184 ^ 

^ ^ 2189 2 166 4 2917 7072 ^ 
P7 = Pj, = 71^ 71^ 1 7i" 

^^ ^^ 1350 675 1080 637875 

641 2 292 4 5119 12448 . 

Pg = 71 71 1 71 

^^ 225 675 1080 637875 
^ ^ 253 2 371 544 ^ 64 , 

p,(, = p,, = 71^ -I 7t"-| 71^ 

Fio Fii ^3^ ^^20 42525 225 

^ ^ ^ ^ 343 751 2 19 4 272 g 



Although at the intermediate stages nontrivial numbers like Euler's y occur, the final 
results are rational polynomials of degree 3 in 71^ and hence are simply expressible by 
use of Riemann's zeta function at even integers 



TT 7C TT 



In the remainder of this section we describe the computational strategy. 

Firstly, we note that not all of coefficients are independent but some of them are 
related, reflecting the diagrammatic symmetry, 

P5(^l,^2,^3)=p4(-^l,-^2,-^3) pll(^l,^2,^3)=Pl0(-^l,-^2,-^3) 

Pl5(^l,^2,^3)=Pl2(-^l,-^2,-^3) pl4(^l,^2,^3)=Pl3(-^l,-^2,-^3) 

P3(^l,^2,^3) = p2fe,^2,^l) P8(^l,^2,^3) = p7fe,^2,^l) 

Pl4(^l,^2,^3)=Pl2fe,^2,^l) Pi5(^l,^2,^3)=Pl3fe,^2,^l). (HO) 



The first four relations hold due to the up-down symmetry and the last four ones 
due to the symmetry w.r.t. reflection at the anti-diagonal line. This reduces our task 
considerably. The explicit expressions of py (pj) in terms of co functions are still very 
much involved. Below we present the simplest case pi for illustration. Supplementary 
arguments are given in Appendix D. The expression for pi contains up to trilinear terms 
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m CO, 



PUSbS2,S3J- ^:^ hCj 0)[C,i ,C,2)+^2 ®lSl)SlJ+C3 (0(.qi ,^2 J 

+ permutations and negation. 

We use the shorthand notation, ^"^ :=^j± i. The last term contains the distinct terms 
under permutations of {^,} and negation of spectral parameters ^,- — )■ — ^;. We regard 
"distinct" after applying the symmetry properties (o(X,-,Xy) = (o(Xy,X;) and (o(X;,Xy) = 

(o(— X,-, — ^j). The coefficients c ' are rational functions of ^1,^2,^3 and their explicit 
forms are listed. 

, , , (^?2+4')(7^i3(^i3-2i)^23fe3-2i) + 80^2^ + 2432) 



cS^ 


(^1 


4' 


(^1 


c^i 


(^1 


cf 


(^1 


4^ 


(^1 


cf 


(^1 


cf 


(^1 


cf 


(^1 


c? 


(^1 


c? 


(^1 


ct 


(^1 



,^2,^3) 



X2X. 



3 



1680(^12 - 2i) (^12 + 2i)^i3 (^13 - 2i)^23 (^23 " 2i) 
-7(-^12^23 +2Y- 7(^13^23 + 2? - %X2^U {^2 + \\, + ^^3 + §) " ^056 

35^2^(^i2-2i)(^i2 + 2i)^2^(^i3-2i)(^i3 + 2i) 

(^12 + 2i)(^12-6i) (7(^13^23 +2)^ + 2^12(71^13^23 +40(^12-41)) +2084) 

-.2 



,^2,^3) 
,^2,^3) 
,^2,^3) 
,^2,^3) 
,^2,^3) 



1680^2^^13(^13-21)^23^3+21) 

(^23 + 4i)(-^23+4i);.S'' (^1,^2,^3) 



12600^2^ (^12 - 21) (^12 + 21)^2^ (^13 - 21) (^13 + 21) (^23 - 21) (^23 + 21) 

/(^12 -41) (^12 + 41) (^23 -41) (-6^2 + 832e2 + ei - 50^2^^2^ - 1520^21^23 - 27328) 

3150^12(^12 - 21) (^12 + 21)^13(^13 - 21) (^13 +21)^2^ 
(^23 +21)4') (^1,^2,^3) 



12600^2^(^12 - 21) (^12 +21)^2^(^13 -21) (^13 +21)^2^ 

(^13 -4i)(^23 + 4i);.f (^1,^2,^3) 

31501^2^(^12 + 21)^2^^2^ 

(^i2-4i)(^i2+4i)(^i3-4i)(^i3+4i);.f (-^2,-^3,-^1) 



31501^12(^12 -2i)(^12 +21)^13(^13 -2i)(^i3 +21)^2^(^23 -21) 

,^2,^3) = ^4'^(^i, ^2,^3) 



,^2,^3) = ^4'^(^i, ^2,^3) 
,^2,^3) = ^4'^(^i, ^2,^3) 
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37800^2^(^12 - 2i) (^12 +21)^2^(^13 -2i) (^13 +2i)^2^ 
4^^(^1,^2,^3) = 



(3),. . ^ /(-6g2+432e2 + ei-5648)(^i2 + 4i)(^i3-4i)(^23+4i) 

■^^'^"^'^ 9450^12^2^^2^ 

(3) . .^_ 4 (I05e^-840e3+ 23632^2 + 420e2ei + 345856^2 -4172^2 + 2891776) 
c, (^1,^2,^3) --4^25^2^^^^^_2.)^^^^_^2i)^2^^^^3_2i)(^^3+2i)^2^(^,3_2i)(^23+2i) 
(3).. . / (-6^2 + 432^2 + ej - 5648) (^12 - 4i) (^12 + 4i) (^13 - 4i) (^23 - 4i) (^23 + 4i) 

""" ^^^'^"'^'^ 9450^12(^12-2i)(^12 + 2i)^2^^23(^23-2i)(^23+2i) 

(3).. . .^^ ;7P(^i,^2,^3)(^23+4i)(-^23+4i) 

5 ^(,l,q2,<,3; 37800^2^(^i2-2i)(^i2+2i)^2^(^13-2i)(^13+2i)(^23-2i)(^23+2i) 

where 62 and ej, are elementary symmetric functions, defined as coeffcients of t^ and 
P in the expansion of (1 +?^i2)(l +?^23)(1 +^^3i)- The nontrivial polynomials in the 
numerators are defined by 

;7P(^i,^2,^3) = 1748992 + 384e^ + 53248^2 + 106496^^3 

+ ^12^13 (-47662^12^13 + 35^^^^3^ _ 196^2^^23 + 12720^12^13-4800^2^ - 136704) 

/jf (^1,^2,^3) = -57835521- 606208^23 + 384e^^23-2304ie| 

+ ^22^2^ (^23 - 6i) (-476^2 + 35^2^^2^ - 196^12^13) - 9728^2^23 + 2324481^2 
- 80^?2^?3 (37^23 +562i) - 320^12^13 (l5^^3 " 901^23 " 76^23 - 2888i) 

/.f (^1,^2,^3) = 22592 - 32801^12 -6ei - 88^2 - 60i^i2 (^?3 +^23) 

+ ^31^32 (l640 + (^li + 101^12 + 30) ^31^32^ 

/jf '(^1,^2, ^3) = 5783552- 384 (e^ + 3e^)+e^(l40^?2^?3 + 33792) +5184^2^12^13^23 -803840^2 
+ ^12^13 (^12^13 (35^2 + 336^^3 + 13304^2^ + 13472) 
+ 70^?2^?3 + 1296^2^^2^ _ 169728^2^ - 1495552) . 

The above result is checked against direct computations (of some density matrix ele- 
ments ) for fixed Trotter numbers. 

7 Conclusion 

The present report offers a novel way to approach the exact evaluation of correlation 
functions of higher spin su{2) chains in a quantitative manner. It utilizes the discrete 
functional relations of qKZ equation type, as well as the direct fusion procedure. 
Explicit results for shorter segments (m = 2,3) clearly show an improvement of the 
present approach compared to the former formulation in [16]. The correlation functions 
are not given by complex contour integrals but they are given in simple factorized 
forms, using rational functions of inhomogeneities and only one nontrivial ingredient CO. 
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One of the direct consequences of such nice expressions is exemplified in Section 5, 
and Section 6: correlation functions of the S = I chain contain only Riemann's zeta 
function with even integer arguments. This complements the famous conjecture in [8] 
that correlation functions of the S = ^ chain are given in terms of zeta functions with 
odd integer arguments. 

Clearly many subjects are still left open. The appearance of the zeta function 
with odd integer arguments for the S = j chain is naturally explained via the explicit 
construction of the solution to the qKZ equation using ^-oscillators. We believe that an 
analogous reasoning is necessary and possible for the higher spin case. 

It becomes clearer that bulk quantities (e.g. the specific heat) of some novel mate- 
rials, e.g., spin-ladders and spin-nanotubes are described by models with higher rank 
Lie symmetry. Their quantum correlations are definitely the next promising goal of our 
understanding of these materials. Therefore the extension of the present approach to 
the su(«) case is highly desirable. 

We hope to come back to these issues in the near future. 
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M. Jimbo, A. Kuniba, T. Miwa for stimulating discussions. AK is grateful to Shizuoka 
University for hospitality. His work was supported by JSPS. JS is supported by a 
Grant-in-Aid for Scientific Research No. 20540370. 

Appendix A: Auxiliary functions for spin 1 

It is sometimes more convenient to deal with polynomials rather than with rational 
functions. For this reason a different normalization of the elementary /^-matrix was 
used in [34]. This leads to differently normalized transfer matrix eigenvalues. In order 
to simplify the comparison with [34] we define the functions 

Ai(X)=^_(X-3i)^+(X + 3i)cp_(X-2i)cp+(X + 2i)-A['l(X), (A.la) 

A2 (X) =^_ (X - 4i)^+ (X + 2i)^_ (X - 2i)^+ (X + 4i) • 

(p_(X-3i)cp+(X + i)(p_(X-i)(p+(X + 3i)-APl(X). (A. lb) 

Then, following [34], we introduce 

Xi(X) = e-T^(X-i)^(X-3i)(p(X)(p(X-2i) ^^,^/^'\ (A.2a) 

q{X-l) 

X3 (X) = eT ^(X + 3i)^(X + i)(p(X + 2i)(p(X) ^%^ . (A.2c) 

^(A + i) 

where 

^(X):=^+(X + i)^_(X-i), (p(X):=(p+(X + i)(p_(X-i). (A.3) 

It follows that 

A2 (X) = Xi (X) + X2 (X) + X3 {X) . (A.4) 
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The basic auxiliary functions for spin 1 are 

X,{X + i)+X^{X + i) X2(X-i)+X3(X-i) 

^^^^ = hix^) ' ^^^^ = Ml^T^ ' ^^-^^ 

with corresponding capital functions 

5S(X) = l + b(X), ^(X) = l + b(X). (A.6) 

In [34] the nonlinear integral equations (48) were derived from a set of functional 
equations satisfied by the functions b, b, !B, 53 together with 

^W = ^,(X-3.).KX+^iML2.).p(X+2i) - i-W-'+^'W- ('^■') 

In appendix B we present a derivation of integral equations by use of algebraic relations 
exposed below. 

b(A) = e rAi(A)— T , .; s — 7-. ^"^^ r, (A.8a) 

^^ ^ >(X + 4i)^(X + 2i)(p(X + 3i)^(X-2i)' ^ ^ 

A2(X) = 53(X-i)X3(X)=^(X + i)Xi(>i), (A.Bc) 

yfl) = Ai(X-i)Ai(X + i) 

^^ ^(X-3i)^(X + 3i)(p(X-2i)(p(X + 2i)' ^^^ 

Note that from (A.7) and (A. 8c) we obtain 

1 , ,a)-i - ^(^) - Ai(X + i)Ai(X-i) 
+^^^^ "K^" Mil) ^^-^^ 



Appendix B: NLIE with straight contour integrations 

To proceed further, it is convenient to consider equations in Fourier space. For a 
smooth function f(k) we define 

Equations (A.8c), (A.2a), (A.2c) and (A.8d), after taking the logarithmic derivative and 
the Fourier transform, read 

dlA2 = e"*^dlB + dlXs = e*^dlB + dlXi , (B.2a) 

dlXi = (&-'' + e-^'')d\<^^+d\(^Q + e-^''d\(^^+e^^d\q+-e-''d\q^, (B.2b) 



dlXs = f e*^ + e^M dl(^+ + dlcp,, + e^Ml(p+ + e^^Mlq^ - eMlq+ , (B .2c) 

dlY = (&'' + &-''] dlAi - e"^*^dl(^_ - e^*^dl(^+ - e^^'^dlcp^ - e2*dl(p+. (B.2d) 
These equations simplify in different ways depending on the sign of k. 
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k>0. 

Here dl(|)+ = dl(p+ = dlq_|_ = 0, hence from (B.2a)-(B.2c) 

dlXs - dlXi = 

Q-kj^Q-3k\ ^y^ _ Q-2k^y^^ ^ (e-'' + e-^^^ dlq_ = -e^'^dlB +e*^dlB, (B.3a) 

from which we can express dlq_ in terms of dl(|)_, dlcp_, dlB and dlB. From (B.2d) we 
obtain dlAi in terms of dl(|)_, dlcp_ and dlY. These expressions are then inserted into 
the Fourier transforms of the logarithmic derivatives of (A.8a) and (A.Sd) yielding 

dlb = dlAo + , ^ , f -e"^dl(b_ + dlY + e"*^dlB - e^^dlB^ , (B.4a) 

dly = dlcpo + -, r f -e"Ml(p_ + dlB + dlB^) . (B.4b) 

e* + e " V / 

The explicit (|) and cp terms yield (note that dl(|)_ refers to the Fourier transform of the 
logarithmic derivative of (|)(A,) in the lower half plane, it is not the Fourier transform of 

e *■ A'^ e"^ — e "'^ 

^^'i^o- k I -k '^^^- = -i^ k , -k ^ ^^-^^^ 

e^ + e" 2e'' + e'' 



dl^o - ^I^dlcp^ = -i^r^e^^^ ( 1 - e'S^) . (B.5b) 



k<0. 

Here dl(|)_ = dl(p_ = dlq_ = 0, and a derivation similar to above yields 



with 



dlb = dl(bo + , ^ , f -e-Ml(|)+ + dlY + e^MlB - cMib) , (B.6a) 

e'' + e " V / 

dly = dlcpo + -^ ^ ('-e"'^dlcp+ + dlB + dlS) , (B.6b) 



g-i \7 ^uk _ p-uk 

^^^0- ^ , _^ dl(|)+ = -i- , _;, , (B.7a) 

gK_l_gK ze^ + e" 



dl^o - ^I^dlcp+ = -i^r^e^^^ ( 1 - e''] . (B.7b) 



These results and an expression for dlb similar to that for dlb are Fourier trans- 
formed, yielding integral equations for the derivatives of the logarithms of the auxiliary 
functions b, b and y in terms of !B, 53 and Y. Integrating these equations and fixing the 
integration constants from the known asymptotics completes our derivation of the NLIE. 

We are interested in Ai for which we derive from (B.2d) 

dlAi = ^ ^ _^ [dlY + e-^l^^ldlc^zp +e-2|^ldl(p^') . (B.8) 
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Appendix C: A simple nearest neighbour correlator 



We are going to relate the derivative of the eigenvalue ^ In AI^I (X) of the spin-1/2 

column-to-column transfer matrix t^^\X;fj.) with spectral parameter X of the spin-1 
system modified by two additional horizontal spin-1/2 lines with spectral parameters 
fj + i and /i — i + 5 to a matrix element of the density operator dI'1 (k, fi) of the unmodified 
spin-1 system for two neighbouring spin-1/2 spaces with spectral parameters X,iu. 

We follow [1]. Be (*P| and |*P) normalized left and right eigenstates of t^^\X;/j) 
with eigenvalue aI^I {x',fj). We then have 
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ln{AW(X;Ai)} 



5=0 38 



ln{(»F|?W(X;/i)|»P)} 



(»F||^W(X;a/)|»P) 



8=0 {^>\tm{x-fi)\^>) 



5=0 

(C.l) 



where for the last equality we used that (*P|*P) = 1 for all 5. Next we use the factoriza- 
tion of (*P| and |*P) at 6 = into eigenstates of the original spin-1 system (without the 
pair of additional horizontal spin-1/2 lines) times a singlet of two spin-1/2 objects. By 
use of the Yang-Baxter algebra and unitarity we find 



o r^ r^ 



U[ 



U\ 
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M/V 



kfl 



UN 



Figure 5: Graphical illustration of the last equation in (C.l). 



35 



In 



aW(X;a/)} 



5=0 



tr{DW(X,A<)/?(X-i,A/-i) |/?(A/ + 5-iA-i)L } 
tr{D[il(X,A/)/?(X-i,A/-i)/?(A'-iA-i)} 



(C.2) 

Using the normalization of the "R matrices used in the definition of f I'l iX',lj) we find 
unitarity in the following form 



R{X-i,^l-l)R{^l + h-i,X-i) 



5=0 



(C.3) 



and the derivative 



R{X-i,ii-i) ^R{^l + b-i,X-l\ 



2i 



5=0 



(X-A/)2 + 4) 



l_p[l] 



(C.4) 
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Upon introducing the function Q.{'k,iLi) by 



2i 



^(^'^) ■= a-j2,4/ '-{^'"(^^A')f'"} , (C.5) 



we find (63) and the equivalent (62). 
Appendix D: More on 3-site density matrix 



In Section 6, the coefficients py were introduced and only the explicit form of pi 
was given. Here we comment on the other ones. 

Although their explicit forms in terms of co are very involved as already commented, 
the expressions in terms of particular 6-site density matrix elements of spin ^ chain are 
still manageable. 

Let us define three elements 



'-'=(«'")::::::. ^.=(o™):::::;. o,=(Dmt:i: 



where the spin-1/2 density matrix D^^' depends on the six arguments {^^ , ^^ , ^ J , ^2 ) ^3^ ) ^3 ) 
where ^^ = ^y it i for j = 1,2, 3. Then we find explicitly 

14(21 + ^1 -^2) „ f t+ .- .- t+ .- . + . 2(-4i + ^2-^3) ^ f .- tt+ tt+ .+ tt- t-, 
9(£ -£2) — ''''^^1 '^^2 '^Sl ,-S2 .-S3 .-S3 )+ 3/ 2i + ?2-£3) *''^^i '^^1 '^'=2 '-S3 '-S2.-S3 , 

on;? ? -LO? ? ? ?2)i X ^ G6(-qi ,-^1 ,-^2 ;-S2 '-S3 '-S3 ) 
"'°'^^^^ + '^'^^^^"^^^^^N(^,-^2)(2l + ^l-^2)(-2l + ^2-^3) 

4(4i + ^l-^2) „ . ,- _.- .+ _.+ ,- + 2(4i + ^2-^3)(-6i + ^2-^3) ^ ..- p+ .- .- .+ .4 
2i + ^l-^2 ^^ ^" ^2' ^" ^2' ^'' ^'^ 3(-2i + ^2-^3)(2i + ^2-^3) ^^"^"^"'^''^2'^' 

3fe-^3)(-2i + ^2-^3) '^^"^"^^'^^'^^'^2^ 3(-2l + ^l-^2)(2i + ^l-^2) '^^"^2'^1'^2'^3'^: 

(-21 + ^1 -^2)(-4i + ^i -^2) ^ ..- .+ .- p+ .- . + . 
3(^l-^2)(2i + ^l-^2) C.6(^2'^2'^l'^l'^3'^3) 

-(48i^i-24^?-96i^2+4^i^2-10i^fe-4^i + 30/^l^i + ^2^i-20i^i-2^1^i + ^^ + 48i^3 
+ 44^1^3 + 101^1^3+4^2^3 - 401^1^2^3 - 2^?^2^3 + 301^1^3 +4^i^i^3 - 2^^^3 - 24^2 + ioi^j^2 

I )i2p2 in;? ?2 ^il p )i2 , ii2ii2-, <J6(qi ,^1 ,C,2 'S2 'S3 'S3 ) 

+ Si S3 - IO1S2S3 - 2SlS2S3 +S2S3)^7^ y w ^- . )i FTTk e \^o; 1 ii S 

9 Si - S2(-2i + qi - ^2(^2 - ^3)(2i + q2 - S. 
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14(2i + ^l-^2) 

9fe-^2) 



G5(-^^-^2,-^^>-^2^-^3">-^l)-|^^|J|^G6(-^^,-^^-^2^-^^-^2,-^3") 



+ 2(56i - 28^1 - 8i^2 _ 4i^j ^^ + ^?^2 + 12i^i - 2^i ^^ ^ ^^ ^ 28^3 + 20i^i ^3 - ^?^3 - 201^2^3 



-2^1^2^3-^i^3 

4(41 + ^1-^2) 
2< + ^l-^2 



G6(-^r.-^^-^2.-^^-^3'.-^3+) 
^(^l-^2)(2i + ^l-^2)(-2i + ^2-^3) 

G(,{-^i,-^2'^^l '^^2 '^^3^ '^^3 ) 



+ (481^1 - 24^? - 96i^2 +4^i^2 - 10i^?^2 - 4^i + 30i^it,l + ^f^^ _ 20i^3 _ 2^j^3 ^^4 

;^^3 



+ 48i^3 + 44^1 ^3 + 10i^?^3 + 4^2^3 - 40i^i ^2^3 - 2^f ^2^3 + 30i^i^3 + 4^i ^i^3 - 2^^^3 



94?2 I in;? ?2 I )i2e2 ,ni??2 ^)l t )i2 , ii2ji2\ ^ <J6(si .Si 1S2 'S2 'S3 -S3 ) 

-24^3 + 101^1^3+^1^3 - IOIS2S3 - 2^i^2S3 +S2S3) X W7i g w ^; , r ^TTp ;i \/o^_Li^ TT 

9(Sl-S2)(-2l + ^l-^2)(S2-S3)(2l + ^2-S3) 

2(41 + ^2 -^3)(-6i + ^2 -^3) „ ..- .+ .- .- .+ ,+ (2i + ^2-^3)(4i + ^2-^3) ^ ,.- .+ .- .+ .- . + . 

+ 3(-2i + ^2-^3)(2i + ^2-^3) '^^"^"^^'^^'^2'^^^+ 3fe-^3)(-2i + ^2-^3) "^^"^"^^'^^'^2'^^^ 

2(-4i + ^i-^2)(6i + ^l-^2)„ ,.+ .+ .- .- .- a+^ , (-21 + ^1 -^2)(-4i + ^1-^2) ^ ,.- .+ .- .+ .- . + . 

><J-6lSi ,S2 .Si >S2.S3 >S3) + TrT. FT^^TT? FT; <^6lS2 .S2 >Si ,Si ,S3 .S3 ) 



3(-2i + ^l-^2)(2i + ^l-^2) 

" "eisi iSi )S2 1S2 'S3 1S3 



3(^1 -^2)(2i + ^1-^2) 



P4(^l,^2,^3) = ^?^(^l,^2,^3) 



14 



G5(-^r.-^2'-^2+.-^3"'-^^)-2^6(-^r'-^^-^2'-^2+.-^^-^3") 



4(-2i^l+4i^2 + ^1^2-^i-2i^3-^1^3+^2^3) 



2(41 + ^2-^3) 
^2-^3 



(^1-^2) (^2 -^3) 
-^6(^1 Xi ,^2 '^3 '^3^'^2 



F6(^+,^r'^2'^^^3^^3+) 



2(-41 + ^l-^2) 
^1-^2 



M^t^^l M^^2 '^3 '^3 



^ 14(21 + ^1 -^2)G6(-^^-^2'-^r'-^2+'-^3''-^3+) I 2(-4/ + ^2- ^3)G6(-^r'-^n-^2+'-^^'-^2'-^3" 



9(^1-^2) ' 3(-21 + ^2-^3) 

-2(561-28^1 -81^2_4i^j^2^^2^^^121^2_2^j^2^^3^28^^^20i^j^^_^2^^ 

on:? ? _L9? ? ? ?2)i \ ^ G(,{-c,y ,-^1 ,-^2 '-S2 '-S3 '-S3 ) 
"'°'^^^^ + '^'^^^^"^^^^^N(^l-^2)(2l + ^l-^2)(-2l + ^2-^3) 
- (481^1 - 24^? - 961^2+4^1^2 - 101^1^2 - 4^1 + 30i^i^^ + -^A - 201^^ - 2^1^^ + ^^ 
+ 481^3 + 44^1 ^3 + 101^?^3 + 4^2^3 - 401^1^2^3 - 2^f ^2^3 + 301^i^3 + 4^i ^^^3 - 2^^^3 - 24^^ 



2s2n 



+ loi^i^^ + i,n - 101^2^^ - 2^i^2^^ + m) X 



Gei^i '^1 '^2 '^2 '^3 '^3 . 



9(^l-^2)(-21 + ^l-^2)(^2-^3)(21 + ^2-^3) 



2(41 + ^2 -^3) (-61 + ^2 -^3; 



G6(^r'^r'^2'^3^^^^i: 



4(41 + ^1-^2) 
2/ + ^1-^2 



G6(-^l ,-^2 '^^1 '^^2 '^^3 '^^3 ) 



3(-21 + ^2-^3)(2i + ^2-^3) 
(21 + ^2-^3)(41 + ^2-^3) ^ ,.- .+ .- .+ .- ,+, 2(-41 + ^i-^2)(61 + ^l-^2) ^ ,.+ .+ .- .- .- . + . 
3(^2-^3)(-21 + ^2-^3) '^^1 '^' '^^ '^^ '^2 '^2 ) 3(-2l + ^i -^2)(2l + ^l -^2) "^^l '^2 'i^l '^2 '^3 '^3 ) 

(~2l + ^l-^2)(-4l + ^l -^2) ^ /J,- tt+ E- j,+ tt- J,+ ^ p .e- (1+ 11- p+ 11- p+N 

3(£ _£2)(21 + £ -to) -^^ ''^2 'Si 'Si 'S3 'S3 )--n5lSi 'Si 'S2'S2'S3 'S3 ) 
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-2F6 



2(4i + ^2-^3)^6(^^^^.^2^^3'.^^^2+) 2(-4i + ^l-^2)F6(^+.^^^^^^2■^3~.^■ 
-2(561-28^1 -8i^?-4i^l^2 + 



-^r, -^1,-^2,-'- ^2, -^^ -^3 



14(2i + ^i-^2)G6(-^^-^2.-^r>-^2+^-^3'>-^3+) 2(-4/ + ; ,, , _, ,. ,^ 
-^"^ ^ ^ 3(-2i + ^2-^3) 



^?^2 + 12i^i-2^1^i + ^^ + 28^3 + 20i^i^3-^?^3 
C-?- -P+ -Pr -P+ -?r -P+^ 



-';2 

^2-^3)G6(-^^:-^^^^2+.-^^-^2.-^3 



on;? ? ^o? ? ? ?2e ■, ^ <J6( Si 1 Si ; S2 > S2 : S3 ; S3 ) 

20i^2^3+2^,^2^3-^2^3) X 9(^,_^^)(2, + ^^_^2)(-2: + ^2-^3; 

+ (481^1 -24^?-96i^2+4^i^2-10i^?^2-4^i + 30i^i^i + ^2^i-20i^3_2^^^3^^4^4gi^^ 

+ 44^1^3 + 10i^?^3 +4^2^3 -401^1^2^3 -2^?^2^3 +301^1^3 + 4^1 ^2^3 -2^i^3 -24^2 

^ifiiP ?2 , )i2)i2 ,n;? ?2 ^t t )i2 , ii2e2-,^ '^eVSi .Si 1S2 >S2 'S3 ;S3 ) 

+ 101^1^3 + ^1^3 - 101^2^3 -2^1^2^3 +^2^3) X 9(^,_^,)(_2i + ^i-^2)fe-^3)(2l + ^2-^3) 

4(41 + ^1 -^2)G6(-^r. -^2. -^^-^2+.-^3''-^^) , 2(4i + ^2-^3)(-6i + ^2-^3) ^ ,.- .+ .- .- .+ ., 
+ 2/ + ^l-^2 +3(-2i + ^2-^3)(2i + ^2-^3)^''^^"^"^2'^^'^2'^^ 

(2i + ^2-^3)(4i + ^2-^3) ^ ..- .+ .- .+ .- p+, , 2(-4i + ^i-^2)(6i + ^i-^2) ^ ,.+ .+ .- .- .- . 
+ 3fe-^3)(-2i + ^2-^3) '^^' '^' '^^ '^^ '^2 '^2 )+ 3(_2i + ^, -^2)(2i + ^i -^2) "^^i '^2 .'^i .^2 .'^s '^: 

^ 3(£i-£2)(2i + ^l-^2) '^<>(^2,^J,^r'^n^3"'^3")+'f6(^r'^n^2'^2"'^3"'^3") 



P7(^l,^2,^3) = ^3^(^l,^2,^3)-^G5(-^r.-^2.-^2+'-^3^-^t 



2(-4i + ^2-^3) p..+ ,- ,- .+ .- + 2(41 + ^2-^3 
1^^7|^^ ^6(^1 ,^1 ,t,2 '^2 .^3 .^3 ) + ^2-^3 

4(21 + ^1-^2) 



■F6(^+,^r,^2,^3-,^+,^2^ 



9(^1-^2) 



G6(-^l ,-^2 '^^1 '^^2 '^^3 '^^3 



2(-4/ + ^2-^3 
3(-21 + ^2-^3: 



■G6(-^^--^^-^2^-^3^-^2>-^3": 



2(-8- 81^1+41^2 + ^1^2-^^41^3 -^1^3 +^2^3) ^ . ,- _,+ .- _,+ ,- .+. 

9(^l-^2)(-21 + ^2-^3) ^'^ ^" ^" ^2' ^2' ^^' ^'^ 

-(481^1 -24^2+481^2 + 76^1^2- 10i^?^2-4^i-61^i^2^^2^2^ 161^3 _2^^^3^^4_9(,i^^ 

- 28^1^3 + 101^?^3 - 68^2^3 + 321^1^2^3 - 2^?^2^3 - 421^2^3 +4^i^^^3 - 2^i^3 +48^? 



26i^i^^ + ^f^^ + 261^2^2_2^^^2^2^^2^2) 



G6(^l ,^1 ,^2 '^2 '^3 '^3 . 



[Gd%-,Xh^2X-,MM)-^^rfi^Y¥^^^i^G,(%-,^i+,k-,Xh^2M) 



9(^l-^2)(-21 + ^l-^2)(^2-^3)(2i + ^2-^3) 
2(41 + ^2 -^3)(-61 + ^2 -^3) 
3(-21 + ^2-^3)(21 + ^2-^3) °^"""""2-3..2>.3/ 3(^2-^3)(-21 + ^2-^3) 

4(-41 + ^i -^2)(6i + ^i -^2)G6(^^^2^^^>^2 .^3'.^^) , 2(-21 + ^i -^2)(-41 + ^i -^2) ^ ,.- ,+ ,- ,+ ,- ,+. 

+ 3(^1 -^2) (21 + ^1-^2) ^6i^2:';2.^1.^1.^3.^3) 



-/'6(^r-^^^2-^^^3^r 



3(-21 + ^i-^2)(21 + ^i-^2) 
3+) 
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P9i^U^2X3) = ^miX2,^3)-lGs{-^i,-^^,-^+,-^^,-^+)-2F,{-^l,-^+,-^^,-^+,-^+,-^-, 



2 (-41^1 +81^2 + ^1^2-^1 -41^3-^1^3 +^2^3 
(^l-^2)fe-^3) 



-F.i^t^^i^^lM^^iM) 



2(4i + ^2-^3)p fp+ t- J,- p- p+ s+\ I 2(-4i + ^i-^2) p /;,+ J,- p+ p- p- t+- 
^b[<ii ,?! ,<;2'';3'S3 'S2)^ tt E: ^6(';2'Si ,<ii ,<i2,^i,^i. 



^2-^3 
4(2i + ^l-^2) 



^1-^2 



9(^1-^2) 

-2i^f + 36! 

■14i^2^3 + 14^1^2^3-7^i^3)x 



G6(-^l ,-^2 '^^1 '^^2 '^^3 '^^: 



2(-4i + ^2-^3 



3^^3(-2i + ^2-^3) 
4^1 ^i + 7^^ + 
G6(-^l ,-^1 ,-^2 '~^2 '~^3 '~^3 ) 



G6(-^r--^^-^2^-^3+--^2>-^3") 



(321-52^1 -2i^K36^2-10i^l^2 + 7^?^2 + 121^1- 14^1^1 + 7^^ + 16^3 + 141^1^3 -7^?^3 



9(^1 -^2) (21 + ^1-^2) (-21 + ^2 -^3) 
-(481^1 -24^?-96i^2+4^i^2-10i^?^2-4^i + 30i^i^2^^2^2_20i^3_2^j^3^^4 

+ 481^3+44^1^3 + 101^2^3 +4^2^3 -40i^i^2^3 -2^fe^3 + 30i^i^3 +4^1 ^i^3 -2^i^3 -24^2 

, in;? s2 , )i2)i2 ,n;? ?2 911 11 )i2 , )i2ii2-,^ ^6(^1 .q, ,^2 .S2 '^3 'S3 ) 

+ 10:^1^3 + ^1^3 - 101^2^3 -2^l^2^3 +^2^3) X 9(^^_^,)(_2: + ^i-^2)fe-^3)(2i + ^2-^3) 

+ 2i + ^l-^2 '^ ^" ^2' ^1' ^^2. ^3. ^3) 3(-2i + ^2-^3)(2i + ^2-^3) ''^^"^"^2'^3'';2.^3) 

(2i + ^2-^3)(4i + ^2-^3) ^ ,.- .+ .- .+ .- K+^ 2(-4i + ^i-^2)(6i + ^i-^2) ^ ,.+ .+ .- .- .- . + . 
3fe-^3)(-2i + ^2-^3) '^^' '^' '^^ '^^ '^2 .^2 ) 3(-2i + ^i -^2)(2l + ^l -^2) "^^l '^2 .^1 .^2 .^3 >^3 ) 

(~2i + ^l — q2)( — 4i + qi —^2) ^ ,r- ji+ e- 11+ t- p+-. „ /e- ji+ i^- p+ r- P+-. 

3(Pl-p2)(2i + £l-b) '^*>(<;2'S2'Si ,«;i ,';3'S3i-'n5(<;i '^i ,<;2'S2'S3.S3J 



Pl0(^l,^2,^3) = ^?^(^l,^2,^3)-^G5(-^r,-^2,-^+.-^3"'-^3+) 

2 (-4/^1 +81^2 + ^1^2- ^1-41^3-^1^3 +^2^3) ^.,+ . . .+ . . + . 

^2(4i+^2_^^^ 2(-^+|p^^^^ 

S2 - S3 SI - S2 

9(l~~^l2)~ '-S2.-S1 .-S2.-S3 '^^3)^ 3(-2i + E2-^3) ' '^S2.-S3 .-S2>-S3) 

2(-8+10l^l-14i^2 + ^1^2-^i + 4i^3-^1^3+^2^3) ^ ,.- .+ .- .+ .- . + . 
9(^,-^2)(-2l + ^2-^3) ^'^ ^" ^" ^2' ^2' ^3. ^3) 

-(48i^i-24^?-96i^2+4^i^2-10i^?^2-4^i + 30i^i^i + ^2^i-20i^3_2^^^3^^4 

+ 481^3 + 44^1 ^3 + 10i^?^3 + 4^2^3 - 40i^i^2^3 - 2^?^2^3 + 30i^i^3 + 4^1 ^i^3 - 2^3^3 - 24^2 

+ 101^1^3 + ^1^3 - 101^2^3 -2^l^2^3 +^2^3) X 9(^^_^,)(_2i + ^,-^2)fe-^3)(2l + ^2-^: 

2(41 + ^2 -^3)(-6i + ^2 -^3) ^ ,a- i+ i- a- i.+ S4 
3r-?.i+eo-£.V?.i+Eo-£,^'^''^'^l'^l''^2'^3'^2'';3 



3(-2i + ^2-^3)(2i + ^2-^3) 



(2i + ^2-^3)(4i + ^2-^3) ^ ,y- p+ p- .+ tt- t + s 
3fe-^3)(-2i + ^2-^3)^''^^"^"^^'^^'^2'^2) 



2(-4i + ^1-^2) (61 + ^1-^2) ^ ,.+ .+ .- .- .- ).+, (-21 + ^1 -^2)(-4i + ^1-^2) ^ ..- .+ .- .+ .- . 
3(-2i + ^i -^2)(2i + ^i -^2) '^^' '^2 .^1 '^2 '^3 .';3 ) 3(^j -^2)(2i + ^i -^2) "^^2 '^2 -'^l '^1 '^3 '^: 

-/'6(^r.^l'^2'^^^3"'^3+) 
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Pl2(^l,^2,^3) = -^?^(^l,^2,^3) + ^G5(-^r,-^2,-^+,-^3-,-^3+) + 2F6(-^7,-^+,-^2,-^+.-^^-^3" 



(^l-^2)fe-^3) F6(^+,^r>^2-^^^3^^3+) 

-'^'^^J^^^6(^^^^-^2-.^3--^3^^2+)-'^=^^t^^ 

4(2i + ^l-^2) „ . .+ t- .- .+ .- . + x 
~9(^~^^2)~ ^^'il'^'il '^S2.-S3 .-S3 J 

-(32i-52^i-2i^K36^2-10i^l^2 + 7^?^2 + 121^1- 14^1^1 + 7^^ + 16^3 + 141^1^3 -7^?^3 

14i? ? _L TzL? ? P 7p2!i \ „ G6(-Sl .-Si 1-S2 .-^2 '~S3 . -S3 ) 

- 14iq2S3 + 14qiq2S3 -7^2S3J X ^Tv e ^r^; 1 S^ p N/ ^; , g ^T 

9(S1 -S2 2l + ^i -^2-21 + ^2-^3 



2(4i + ^l-^2) 



G(,{-i,i ,-^2 '^^1 '^^2 '^^3 '^^3 



2' + ^l-^2 ""' ^" ^2' ^" ^2' ^^' ^^' 
+ 2(24i^i - 12^2-481^2- 16^1^2- 14i^?^2 + 16^i + 24i^i^2^5^2^1- lOi^i- 10^1^3 ^5^1 

+ 241^3+40^1^3 + 14!^?^3 - 16^2^3-201^1^2^3 - 10^?^2^3 + 61^1^3 +20^1 ^1^3 - 10^1^3 - 12^? 



c2si2 , ^;ji ji2 ,aS S s2 , j:Si2ji2\ ^6(^1 ,C,^ ,^2 'S2 »S3 1S3 ) 

'" 9(^l-^2)(-2i + ^l-^2)fe-^3)(21 + ^2-^3) 



- 41^1^2 ^ 5^2^2 ^41^2^2 _ 10^^ ^^^2 ^ g^2^2 



2(-4.- + ^2-^3)G6(-^r>-^^-^^-^^-^2.-^3') , 8(-3l + ^2-^3)(4i + ^2-^3) ^ ,.- .+ .- .- .+ . + . 
3(-21 + ^2-^3) +3(-2i + ^2-^3)(2i + 52-^3) '^^"^"^^'^^'^2'^^^ 

(2; + ^2-^3)(4i + ^2-^3) ^ ..- .+ .- .+ .- p+^ , 2(-41 + ^i -^2)(6i + ^i -^2) ^ ,.+ .+ .- .- .- . + ^ 
+ 3fe-^3)(-2i + ^2-^3) '^^' '^' '^^ '^^ '^2 '^2 ^+ 3(-2l + ^i -^2)(21 + ^l -^2) "^^l '^2 -'^l .^2 >';3 '^3 ) 

H 3(£ -£2)(21 + 4 -^2) ^<>(^2,^J,^r,^i^,^3",^;[)+-P6(^r,^n^2 ,^^,^3",^;^) 

The rest is easily reproduced by applying the symmetry relations (110). The explicit 
forms of Fa,Ga and Pa in terms of (o(^;, ^y) are too lengthly to be reproduced here. 
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